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Abstract 

Disease surveillance in wildlife populations presents a logistical challenge, yet is critical in 
gaining a deeper understanding of the presence and impact of wildlife pathogens. Erinaceus 
coronavirus (EriCoV), a clade C Betacoronavirus, was first described in Western European 
hedgehogs ( Erinaceus europaeus) in Germany. Here, our objective was to determine whether 
EriCoV is present, and if it is associated with disease, in Great Britain (GB). An EriCoV-spe- 
cific BRYT-Green® real-time reverse transcription PCR assay was used to test 351 samples of 
faeces or distal large intestinal tract contents collected from casualty or dead hedgehogs from a 
wide area across GB. Viral RNA was detected in 10.8% (38) samples; however, the virus was 
not detected in any of the 61 samples tested from Scotland. The full genome sequence of the 
British EriCoV strain was determined using next generation sequencing; it shared 94% iden¬ 
tity with a German EriCoV sequence. Multivariate statistical models using hedgehog case 
history data, faecal specimen descriptions and post-mortem examination findings found no 
significant associations indicative of disease associated with EriCoV in hedgehogs. These find¬ 
ings indicate that the Western European hedgehog is a reservoir host of EriCoV in the absence 
of apparent disease. 


Introduction 

The Betacoronavirus genus includes multiple bat-associated viruses and several human patho¬ 
gens, including severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle 
Eastern respiratory syndrome coronavirus (MERS-CoV) [1, 2], In 2013, a novel clade C 
Betacoronavirus was discovered in the Western European hedgehog ( Erinaceus europaeus) 
in Germany [3]. The virus was subsequently detected in the same species in France [4]. 
Characterisation of these Erinaceus coronavirus (EriCoV) nucleotide sequences revealed 
high nucleotide identity to MERS-CoV [3], the cause of an acute respiratory syndrome in 
humans with high case fatality rates [5, 6]. 

The structural spike protein is key to coronavirus virulence and host specificity [7] and the 
gene that encodes it plays a pivotal role in coronavirus diversity [8]. Analysis of the receptor¬ 
binding region in the S protein of EriCoV revealed that it shared only 36.7% amino acid iden¬ 
tity with that of MERS-CoV [3]. However, there have been no further published EriCoV whole 
genome sequencing attempts, and many questions remain over the origins and evolution of 
EriCoV [2]. Whilst there is no known threat to public health from EriCoV, further investiga¬ 
tion of the virus is warranted to gain a fuller understanding of the impact of coronavirus infec¬ 
tion on wildlife populations. 

Bats are widely recognised and accepted as the natural reservoir hosts of a variety of 
coronavirus species, including SARS-CoV [9, 10]. Both bats and camels have been cited as pos¬ 
sible reservoirs of MERS-CoV [6, 11]. It is considered that bats are the likely natural reservoir 
of MERS-CoV [2, 5], and that dromedaries ( Camelus dromedarius), and potentially other spe¬ 
cies, act as an intermediate host for human transmission [6, 12, 13]. MERS-CoV-like virus 
sequences have only been detected in one other species, E. europaeus [3], Hedgehogs are mem¬ 
bers of the order Eulipotyphla, which is closely related to Chiroptera, the order of bats [14]. 
This close association is what initially led to hedgehogs being investigated as a potential 
coronavirus reservoir [3]. 

Many animal species seem to have the capacity for coronavirus infection in the absence of 
apparent disease, including bats [15], aquatic birds [16] and rabbits when inoculated with 
MERS-CoV [17]. However, dromedaries infected with MERS-CoV develop pyrexia with 
mild upper respiratory signs, consisting of serous to purulent nasal discharge [18]. Lesions 
in infected dromedaries are comparable to those seen in the common cold in humans, includ¬ 
ing inflammation and loss of pseudostratified epithelial cells in the upper respiratory tract. 
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Coronaviruses that affect the gastrointestinal system, such as por¬ 
cine transmissible gastroenteritis virus, typically cause similarly 
mild lesions, including destruction of enterocytes and villus atro¬ 
phy [19]. 

Although still widespread throughout much of Western and 
Northern Europe, the Western European hedgehog population 
in Great Britain (GB) has seen a rapid decline in recent decades, 
with an estimated population decrease of 95% between the 1950s 
and 1990s [20, 21]. The reasons for this decline are not fully 
understood, and it is likely to be multi-factorial (e.g. habitat loss 
and fragmentation, road traffic death and predation) [21, 22]. 
The possible role of infectious disease also needs to be investi¬ 
gated, including the consequences of infection with newly- 
discovered viral agents. Here, we identify the presence of 
EriCoV in British hedgehogs, sequence the full genome, explore 
the spatial distribution of infection and investigate if there is 
any association of EriCoV with disease. 

Methods 

Sample collection and preparation 

A total of 351 samples from individual hedgehogs were tested in 
this study. Of these, 225 were voided faecal samples collected 
from live hedgehogs that had been found injured, sick or aban¬ 
doned juveniles and were consequently admitted to wildlife centres 
for treatment and rehabilitation in 2014 and 2015. A total of nine 
wildlife centres across GB collected and supplied samples, using the 
methods described by Sangster et al. [23] . A fresh faecal sample was 
collected non-invasively shortly after admission and associated 
metadata, including clinical signs, signalment, location and history, 
were collected where available. Faecal samples were stored at 4 
or —20 °C for a maximum of 1 week at rehabilitation centres 
prior to submission. A subsample of each sample collected by 
Sangster et al. [23] was used in the current study. An additional 
four voided faecal samples were collected from the wild as part 
of a Greater London-based hedgehog population survey conducted 
in 2014 to give a total of 229 faecal samples from live animals. 

The remaining 122 samples were of distal large intestinal tract 
contents collected from dead hedgehogs during post-mortem 
examination (PME) from 2011 to 2015 and comprised animals 
with infectious and non-infectious (e.g. trauma) causes of death. 
These hedgehogs had either died in wildlife centres within 48 h 
of admission or had been found dead by members of the public. 
Where possible, PMEs were performed on fresh carcases, but fro¬ 
zen carcases that had been archived at —20 °C were also examined. 
Systematic examination of external and internal body systems was 
performed using a standardised protocol, followed by microbio¬ 
logical, parasitological and histological investigations where indi¬ 
cated by macroscopic findings and when permitted by the state of 
carcase preservation [24]. The contents of the distal large intes¬ 
tinal tract were sampled to most closely reflect freshly voided 
faeces, as the latter had contained the highest EriCoV viral titre 
in a previous study [3]. 

Clinical data and case history 

Faecal sample submission forms provided information on: admit¬ 
tance date, sample collection date, the location where the hedge¬ 
hog was found, sex, live body weight, faecal characteristics 
(including colour and consistency) and reason for admittance. 
PME records were available for 121 of 122 hedgehogs from 
which samples of distal large intestinal contents were tested. 


These comprised data on the location where the dead hedgehog 
was found, the date the carcase was found, sex, age class, carcase 
weight, subjective assessment of body condition (based on muscle 
mass and body fat deposits), gross PME findings and the results 
of ancillary diagnostic tests. 

Viral RNA extraction 

RNA stabilising agent RNAlafer (Ambion, Life Technologies 
Europe, NL) was added to, and mixed with, a pea-sized amount 
of each sample (faeces or distal large intestinal tract contents) 
in a 1.5 ml microcentrifuge tube. Each sample was then stored 
at either —20 or —80 °C. Faeces were homogenised and pooled 
(five samples per pool). Each 140 pi of total pooled sample was 
added to 560 pi buffer AVL-carrier RNA solution and extracted 
immediately following the spin protocol alongside negative 
extraction controls. A guanidine and column-based method 
(QIAamp Viral Mini Kit Qiagen, Hilden, DE) was used for 
RNA extraction, according to the manufacturer’s protocol. RNA 
was then eluted in 60 pi of buffer AVE and stored at —80 °C. 

Detection of EriCoV RNA 

A specific oligonucleotide primer set (EriCoVR and EriCoVF, see 
Supplementary Table SI) was designed based on the EriCoV 
RNA-dependent RNA-polymerase (RdRp) gene sequence of the 
German virus sequence (KC545383) with OligoArchitect™ 
(Sigma-Aldrich*, St Louis, MO, USA) and used with the 
BRYT-Green based GoTaq* 1-Step real-time reverse transcription 
PCR (RT-PCR) system kit (Promega Corporation, Madison, WI, 
USA) to detect EriCoV RNA in an Applied Biosystems 7500 Fast 
Real-Time system. A 10~ 2 dilution of extracted RNA from a posi¬ 
tive EriCoV sample (R618/14, previously confirmed through 
sequencing) was used as the positive control. Nuclease-free 
water (Promega Corporation) was used as the negative control. 
Following the amplification stages, melt curve analysis was per¬ 
formed in order to verify the specificity of the amplicons [25], 

In addition, a 110 bp EriCoV ultramer oligonucleotide was 
used to quantify the limits of detection and assess the sensitivity 
of the designed assay. This was a sequence within the EriCoV pri¬ 
mer RdRp target region, but with the removal of 10 nucleotides at 
the centre to differentiate it from wild-type EriCoV (see 
Supplementary Table SI). A 14-fold dilution series was performed 
and a standard curve was generated using GraphPad Prism 6 
(GraphPad Software, CA, USA). To assess assay specificity, 
RNA of the gammacoronavirus, infectious bronchitis virus 
(IBV), was also included as a control. 

To assess the primer dynamic range and potential RNA deg¬ 
radation, a dilution series of RNA extracted from faeces was 
tested. A ‘spiked’ dilution series was set containing 25% (5 (il in 
20 pi) faecal RNA extraction from an EriCoV-negative pool and 
a subset of negative samples (« = 11) was tested using primers 
for the 18S ribosomal RNA subunit. All samples extracted from 
faeces produced positive results confirming the presence of ampli- 
fiable RNA (data not shown). 

Real-time RT-PCR assays using the EriCoV primers were car¬ 
ried out for all pooled samples. Single samples from positive pools 
were then extracted and tested individually. A subset of amplicons 
from individual samples (n = 10) was also assessed by gel electro¬ 
phoresis and compared to the positive control. Sanger sequencing 
was conducted on three amplicons with EriCoVR and EriCoVF 
primers using standard protocols. Amplicon sequences were 
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manually checked then aligned against the German EriCoV RdRp 
gene fragment sequence retrieved from GenBank (accession num¬ 
ber KC545385) to confirm their identity (MEGA 6.0). 

Whole genome sequencing 

Whole genome sequencing was performed on RNA extracted from 
one faecal sample collected in 2014 (R618/14) which was identified 
as EriCoV-positive by real-time RT-PCR. Sequencing was per¬ 
formed using synthesis technology (Illumina, San Diego, CA, 
USA). RNA was extracted as before and DNA was depleted from 
this sample using the on-column DNase treatment (RNeasy* plus 
mini kit) following the manufacturer’s instructions (Qiagen); the 
resultant RNA was then eluted in 30 pi molecular grade water 
[26]. cDNA was synthesised from 50 ng RNA using a random 
cDNA synthesis system (cDNA Synthesis Kit, Merck, De), according 
to the manufacturers’ instructions. The cDNA was purified using 
Ampure XP magnetic beads (Beckman Coulter, High Wycombe, 
UK), 1 ng of which was processed for sequencing using the 
Nextera XT DNA sample preparation kit (Illumina). A sequencing 
library was prepared according to the manufacturers’ instructions 
and sequenced using an Illumina MiSeq with 2 x 150 bp paired 
end reads following standard Illumina protocols. 

For confirmation of the insertions and deletions of the result¬ 
ing sequence, seven extra sets of primers were designed (available 
on request). The sequences of these short amplicons were com¬ 
pared with the sequence derived by Illumina sequencing. The 
total reads (12 908 682) were mapped to a reference EriCoV 
sequence from Germany (GenBank accession number 
KC545383). A modified SAMtools/vcfutils script was used to gen¬ 
erate an intermediate consensus sequence in which any indels 
relative to the original reference sequence were appropriately 
called. The intermediate consensus was used as the reference for 
subsequent iteration of mapping and consensus calling. The 
total number of assembled viral reads was 1 217 783 (9.43% of 
the total reads). Despite the low proportion of viral sequence 
detected within the total data set, coverage of the entire genome 
was obtained (average coverage depth of 4546 reads). 

Post-mortem analysis and histopathological investigation 

Gross PME findings were categorised into whether respiratory 
abnormalities (including, but not limited to, pneumonia, consoli¬ 
dation, haemorrhage, congestion, granuloma and high Crenosoma 
sp. burden) and/or gastrointestinal abnormalities (including, but 
not limited to, haemorrhage, inflammation, intussusception, rup¬ 
ture, abscess and faeces of abnormal colour or consistency) were 
detected. Gross PME findings were available for both organ sys¬ 
tems from 94 of the 121 hedgehog examined post mortem and 
were included in these analyses. 

A subset of hedgehogs (N = 9) examined postmortem that 
were EriCoV-positive was selected for histological examination. 
Selection criteria focused on carcases likely to provide the optimal 
state of tissue preservation from those available, prioritising tis¬ 
sues from carcases examined fresh or when frozen, from carcases 
with the least evidence of autolysis. A range of available tissue 
samples from each selected carcase was fixed in neutral-buffered 
10% formalin and was prepared for histopathological examination 
and stained with haematoxylin and eosin using standard techni¬ 
ques. Particular attention was paid to examination of the respira¬ 
tory and alimentary tracts, when available, as these are the organ 
systems frequently targeted by mammalian Co Vs. 


Statistical analysis 

Associations between EriCoV-positive status, case history signal - 
ment and pathological findings were explored. Variables were 
first re-categorised and coded using Microsoft Excel before 
being exported to IBM* SPSS* version 22 for Windows (IBM*, 
Chicago, IL, USA) for statistical analyses. Cross-tabulation and 
Pearson’s j 2 tests were undertaken for the binary dependent vari¬ 
able (EriCoV presence or absence) against all other categorical 
independent variables. In addition, where sample numbers were 
sufficient, binary logistic regression (BLR) was performed to 
include consideration of continuous independent variables, such 
as age, and to quantify the strength of association for multiple 
variables simultaneously (year, region). A P-value of <0.05 was 
considered statistically significant. ArcGIS Pro version 11 software 
(Esri, Redlands, CA, USA) was used to map the distribution of 
hedgehogs for which location data were available. 

Results 

Detection of EriCoV 

Of 351 samples examined, 38 (24 faecal samples and 14 distal 
large intestinal tract samples) were positive for EriCoV, giving a 
positive percentage of 10.8%. The viral copies detected in positive 
faecal samples ranged from 1.15 x 10 8 to 2.7X10 4 per 4 pi 
extracted RNA. The limit of detection for the BRYT* Green I 
based EriCoV primer RT qPCR assay was calculated as 2.68 x 
10 3 copies per 1 pi of extracted RNA sample. IBV RNA at 10 -2 
dilution was not detected by the EriCoV-speciflc assay. 


Whole genome sequencing 

The complete genome of the hedgehog EriCoV from sample 
(R618/14, Genbank accession number: MK679660) is 30173 
nucleotides (nt) in length, which is comparable to those 
sequenced from Germany: KC545383 and KC545386 (30 148 
and 30 175 nt, respectively). Complete genome comparison with 
open reading frames (ORFs) from the German EriCoV 
(KC545383) gave identities ranging from 79% (ORF3b) to 98% 
(ORF6, encoding the envelope protein). The genes encoding the 
spike protein of the British and German EriCoV viruses shared 
93% identity (Table 1). 


Geospatial analysis 

The locations for all hedgehogs with available data (n = 344) are 
displayed in Figure 1, along with their EriCoV status. 
Hedgehogs sampled for the study covered a wide area across 
GB, from Devon in South West England to Northern Scotland. 
Hedgehogs positive for EriCoV were found in Southern 
England, East of England, North East England and Wales. We 
found no evidence of EriCoV in hedgehogs from Scotland, even 
though 17.4% (n = 61) of sampled hedgehogs were submitted 
from this country. There was a bias of submissions from the 
East of England, accounting for 44.4% (« = 123) of the total, 
with Norfolk alone contributing 23.1% (n = 64) of the hedgehog 
samples examined, which is explained by a single, large participat¬ 
ing wildlife centre in this region. Although all regions of GB were 
represented, some appeared to be under-represented; for example, 
only 1.1% (n = 3) of submissions came from North East England. 
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Table 1 . Comparison of open reading frame (ORF) lengths of Erinaceus 
coronavirus from Great Britain (GB) and German hedgehogs (Erinaceus 
europaeus) 


Open reading 
frame (ORF) 

GB EriCov 
protein 
length (AA) 

% AA identity with EriCoV/ 
2012-174/Ger/2012 (GenBank Acc. 

No. KC545383) 

ORFlab 

7156 

95.6 

ORFla 

4461 

94.1 

S 

1328 

92.7 

ORF3a 

104 

89.5 

ORF3b 

140 

79.3 

ORF4a 

81 

81.5 

ORF4b 

225 

83.3 

ORF5 

233 

87.9 

ORF6 (E) 

82 

97.6 

ORF7 (M) 

218 

96.8 

ORF8a (N) 

425 

97.2 

ORF8b 

200 

81.4 


AA, amino acid 


Clinical and pathological associations with EriCoV 
detection 

The majority of hedgehog samples were submitted in 2014 and 
2015. The proportion positive in 2015 was higher (28/144, 19%) 
than in 2014 (10/182, 5%) but there were too few submissions 
from other years to determine if there was an association with 
year of sampling. There was a higher number of juvenile hedge¬ 
hogs included in the study in 2015 (n = 68) than 2014 (« = 24), 
but the proportion of juveniles positive in 2015 (19.1%) was com¬ 
parable to the proportion of adults positive (20.3%). 

As coronaviruses are associated with diarrhoea in other spe¬ 
cies, we tested for the possibility of an association between 
EriCoV status and abnormal faeces. The proportion of samples 
from hedgehogs with reportedly normal faecal colour or con¬ 
sistency, which tested positive for EriCoV (14/102), was not 
significantly different from the proportion of all faecal samples 
with abnormal consistency and/or abnormal colour (8/115) 
(P = 0.099). However, a significant association was observed 
between EriCoV-positive hedgehogs and those with green faeces 
(P = 0.004). A significant association was also observed between 
EriCoV-positive hedgehogs and those with yellow faeces (P = 0.034). 
Colour variation within a single sample of hedgehog faeces was 
not uncommon: 42.6% (n = 20) of hedgehogs with green coloured 
faeces were also reported to have faeces of brown, yellow or black 
colouration. No other statistically significant association was shown 
across the other independent variables, including body condition 
(Table 2). A BLR model testing all but the non-mutually exclusive 
variables also showed no statistical significance. 

The proportion of hedgehogs with respiratory abnormalities 
reported at PME that tested positive for Eri-CoV (5/55) was not 
significantly different from the proportion of those without 
respiratory abnormalities that tested positive (7/39) (P = 0.205). 
The proportion of hedgehogs with gastrointestinal tract (GIT) 
abnormalities that tested positive (5/55) was not significantly dif¬ 
ferent from those that had no GIT abnormalities (6/39) (P = 0.35). 


A statistically significant association was shown between 
EriCoV status and geographical region of origin (P = 0.008) 
(Table 2). The highest proportion of EriCoV-positive hedgehog 
samples were submitted from the South of England (34/217, 
16%); however, BLR showed no significant association (P = 
0.678) between EriCoV infection status and wider region when 
other factors including age and year were included. 

Histopathology 

Histopathological examination was conducted on available tissues 
from nine of the EriCoV-positive hedgehogs in this study that 
were examined post-mortem: results are presented with the signal- 
ment, macroscopic findings and other ancillary diagnostic tests for 
interpretation in Supplementary Table S2. Interpretation was com¬ 
plicated by the state of tissue preservation in the majority of cases, 
limiting the ability to appraise evidence of either a host inflamma¬ 
tory response or the presence of mild abnormalities (e.g. loss of 
pseudostratified epithelial cells in the respiratory tract or intestinal 
villus atrophy, as has been reported with some coronavirus infec¬ 
tions in other species [27]). Examination of the respiratory and 
gastro-intestinal tract identified frequent metazoan parasite infec¬ 
tions of variable severity, sometimes in combination with broncho¬ 
pneumonia. A range of other significant findings was detected, 
supported by ancillary testing, including trauma and yersiniosis. 

Discussion 

We detected EriCoV, a novel clade C Betacoronavirus, in wild 
Western European hedgehogs in GB. This virus had previously 
only been reported from Germany and France; in both instances 
from this same species. Our results, therefore, indicate that 
EriCoV is widespread in this European wild mammal. Full gen¬ 
ome sequence of EriCoV from GB showed 94% identity with a 
EriCoV previously detected in Germany, and 78% identity with 
the virus responsible for MERS [3]. 

There was no evidence to indicate that EriCoV was causing 
clinical disease in hedgehogs. The lack of association with detec¬ 
tion of EriCoV and abnormal faecal consistency, poor body con¬ 
dition (i.e. thin or emaciated), respiratory and digestive tract 
abnormalities suggests EriCoV is unlikely to be a primary patho¬ 
gen in this species. This concurs with studies reporting natural 
CoV infections in wild animal species, including bats [15] and 
aquatic birds [16], which do not manifest as clinical disease. 
This is also a common feature of a reservoir or maintenance 
host [28, 29]. Coronavirus infections of man and other animals 
are known to predispose to secondary infections, for example, a 
study recently isolated feline coronavirus and several other enter- 
opathogens from cats with diarrhoea [30], thus making any asso¬ 
ciation difficult to detect should it occur with EriCoV infection in 
the hedgehog. In the current study, histopathological examination 
was complicated by concurrent parasitic and bacterial infections in 
several cases and limited by freeze-thaw artefact and the state of 
tissue preservation, therefore microscopic lesions caused by viral 
infection could have been missed or obscured. Experimental 
infection studies and the development of in situ hybridisation to 
localise EriCoV in tissues would be worthwhile in the future to 
further elucidate the clinical significance of EriCoV infection. 

Despite the lack of association between EriCoV and faecal 
abnormalities overall, there was an apparent association between 
EriCoV status and green or yellow faeces. Green faeces is not 
an unusual finding in hedgehogs, and can often be associated 
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Fig. 1. Geographic origin of hedgehog samples and EriCoV testing result (ArcGIS). 


with Crenosoma sp., Isospora or Salmonella sp. infection [31]. 
Many of the hedgehogs reported to have green faeces also had fae¬ 
ces of another colour, suggesting that faecal colour change is not 


unusual. The colour change or the reporting thereof may not be 
consistent and therefore this result should be interpreted with 
caution. However, the possibility of a causative relationship 
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Table 2. Associations with EriCoV infection 


Independent variable 

Sub-category Number positive/number tested 

Percentage positive (%) 

P value 

Region 

South England 

34/217 

16 

0.008 

Midlands England 

1/16 

6 


North England 

1/25 

4 


Wales 

2/19 

11 


Scotland 

0/61 

0 


Sex 

Male 

3/40 

8 

0.674 


Female 

16/146 

11 


Age class 

Juvenile 

26/97 

27 

0.035 


Adult 

22/234 

9 


Year 

2011 

0/1 

0 

0.001 


2012 

0/13 

0 



2013 

0/9 

0 



2014 

10/182 

6 



2015 

28/144 

19 


Month 

January 

0/1 

0 

0.219 


March 

0/1 

0 



April 

0/29 

0 



May 

5/68 

7 



June 

6/81 

7 



July 

10/61 

16 



August 

4/32 

13 



September 

5/34 

15 



October 

3/17 

18 



November 

4/15 

27 



December 

1/5 

20 


Faecal sample data only 

Consistency of faeces 

Abnormal (overall) 

10/79 

13 

0.361 


Solid 

14/151 

9 

0.499 


Runny 

3/28 

11 

0.921 


Mucoid 

8/57 

14 

0.263 

Colour of faeces 

Abnormal (overall) 

11/73 

15 

0.090 


Brown 

15/161 

9 

0.470 


Green 

10/47 

21 

0.004 


Yellow 

3/10 

30 

0.034 


Black 

1/25 

4 

0.277 

Abnormal faeces overall 

- 

14/102 

14 

0.099 

Post-mortem examination data only 

Carcase state 

Fresh 

4/53 

8 

0.119 


Frozen 

9/51 

18 


Body condition score 

Emaciated 

3/22 

14 

0.231 


Underweight/thin 

0/23 

0 



Healthy/normal 

7/43 

16 



Overweight/fat 

2/11 

18 



( Continued ) 
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Table 2. ( Continued.) 


Independent variable 

Sub-category 

Number positive/number tested 

Percentage positive (%) 

P value 

Gastrointestinal tract abnormalities 

None 

6/39 

15 

0.350 


Detected 

5/55 

9 


Respiratory tract abnormalities 

None 

7/39 

18 

0.205 


Detected 

5/55 

9 



A P-value of <0.05 was considered statistically significant. 


between EriCoV infection and abnormal faeces colour cannot be 
excluded. 

Significantly stronger amplification of CoV RNA has previ¬ 
ously been described in juvenile and lactating female Myotis sp. 
bats [32], and younger age has been shown to be associated 
with increased coronavirus shedding across various bat species 
[33, 34], Although no evidence for such an association for 
EriCoV in the hedgehog was found in the current study, juveniles 
accounted for only 27.6% of the hedgehogs sampled. Examination 
of a larger sample size of juvenile animals would be worthwhile to 
further explore the possibility of age-structured prevalence. 

Hedgehogs positive for EriCoV were detected over a wide 
region of England and Wales, but not from Scotland even though 
17.4% of all samples tested were submitted from this country 
(Fig. 1). Inferring population prevalence from these data is 
inappropriate due to the sampling strategy, relying on conveni¬ 
ence sampling and the catchment area of rehabilitation centres. 
There is apparent sampling bias since 47% of all hedgehog samples 
were from the East of England. Demonstrating regional differences 
in population prevalence and absence of infection would require a 
much larger, random sample of the population, and information 
regarding the population structure of hedgehogs in GB. 

Annual variation in the proportion of EriCoV-positive samples 
was found. The proportion of hedgehogs sampled in 2015 that 
tested positive for EriCoV was higher than that in 2014; a differ¬ 
ence that was consistent irrespective of sample type or hedgehog 
age class (Table 2). This temporal variability in detection rate is 
consistent with previous multi-year studies in bats where signifi¬ 
cant annual variation in coronavirus detection rate was reported 
[34]. Seasonal variation in contact rates, environmental factors 
and fluctuating levels of CoV immunity in the population have 
all been postulated as causes for such fluctuations [34]. Proving 
viral persistence and establishing critical community size is chal¬ 
lenging for wildlife populations [29] . Whilst there is now a wealth 
of studies investigating bat Co Vs, research is lacking for other 
wildlife CoVs. There is evidence that the gregarious roosting 
habits of some species of bat facilitate CoV transmission [35], 
but the transmission dynamics of EriCoV in the more-solitary 
hedgehog [20] are likely to be very different. Establishing whether 
infection is endemic in British hedgehogs will require longitudinal 
sampling, augmented by age-specific serology [36]. 

Overall, the percentage of EriCoV-positive hedgehogs 
detected, at 10.8%, was lower than that reported previously, 
where 58.9% (148/248) and 50% (37/74) of specimens tested posi¬ 
tive for EriCoV in Germany and France, respectively [3, 4]. This 
difference could be due to assay sensitivity or sample degradation. 
Although the limit of detection of the previously used assays was 
not given, average copy numbers in the faecal samples were well 
over the sensitivity threshold of the assay used here. All three 
studies included hedgehogs in wildlife centres. In the German 
study, no details were given for the timing of sampling. In the 


current and French studies, however, faecal samples were taken 
soon after admission (typically within 48 h) in order to reduce 
the possibility of nosocomial transmission. Dog kennels and shel¬ 
ters are recognised as important environments for maintaining 
enteric canine CoVs, with increased length of stay associated 
with increased probability of CoV infection [37]. However, noso¬ 
comial transmission seems to be an unlikely cause of the differ¬ 
ences in sample percentage of positive samples in this study. 
Another potential factor is that the relatively high percentage of 
positive samples in the French and German studies could be sim¬ 
ply due to a much higher prevalence of EriCoV in the mainland 
Europe hedgehog population: the recruited wildlife centres (from 
which the hedgehogs were sampled) were all located in North and 
North Western France, with one located close to the border with 
Germany [4], More studies are required to determine this. 

Hedgehogs in Northern Europe hibernate for 2-5 months of 
the year [20]. EriCoV was detected in both 2014 and 2015, sug¬ 
gesting that the virus overwinters in the hedgehog population, 
but the method of viral persistence is not known. In the current 
study, no association was detected with the month of sampling 
but to explore this further, longitudinal sampling of individual 
hedgehogs could be conducted, including from animals before 
and after hibernation, to determine if there is persistence of infec¬ 
tion in individuals. 

Deriving full genome sequence is an important step in charac¬ 
terising this virus to assess relationships with other CoVs, includ¬ 
ing pathogens of domestic animals and man. Overall identity to 
the previous coronavirus detected in hedgehogs from mainland 
Europe was high (94%), but lower than the identity between 
two German EriCoV isolates (97%) [3]. Previous studies have 
demonstrated that EriCoV shares a common ancestor with 
MERS-CoV, and bat coronaviruses HKU5 and HKU4, raising 
the possibility of a lineage that has been positively selected for 
and which has evolved with the potential for cross-species trans¬ 
mission. Host tropism in the coronaviruses is largely determined 
by the spike protein and its ability to bind to host cell receptors. 
Both the German and the GB EriCoVs have <40% amino acid 
similarity with MERS-CoV at the crucial receptor-binding 
domain on the spike protein, suggesting that binding to the 
human receptor (dipeptidyl peptidase 4) is unlikely. Virus culture 
was not attempted in this study, but was previously unsuccessful 
and this might be due to the cell tropism of the virus [3]. 

The detection of EriCoV in British hedgehogs is a notable 
finding, and the first report outside mainland Europe. Wildlife 
disease research in GB has greatly increased in recent years, par¬ 
ticularly for popular wildlife species such as the hedgehog, where 
citizen science offers a cost-effective approach to achieving 
national disease surveillance [38]. It is important that this con¬ 
tinues if we are to gain a deeper understanding of endemic and 
emerging pathogens of wildlife and their impact on wildlife, 
livestock and human populations. 
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